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ABSTRACT 

It is estimated that a star brighter than visual magnitude 17 is undergoing a 
detectable gravitational microlensing event, somewhere on the sky, at any given 
time. It is assumed that both lenses and sources are normal stars drawn from 
a standard Bahcall-Soneira model of our Galaxy. Furthermore, over the time 
scale of a year, a star 15th magnitude or brighter should undergo a detectable 
gravitational lens amplification. Detecting and studying the microlensing event 
rate among the brightest 10 8 stars could yield a better understanding of Galactic 
stellar and dark matter distributions. Diligent tracking of bright microlensing 
events with even small telescopes might detect planets orbiting these stellar 
lenses. 

Subject headings: gravitational lensing - dark matter - stars 
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1. Introduction 

The ability of one star to gravitationally amplify the light from another was considered 
previously by Einstein (1936), and developed more generally by, for example, Refsdal (1964). 
Realistic rate estimates for gravitational microlensing by massive parallel photometry were 
made by Paczynski (1986) for the LMC and again by Paczynski (1991) for the Galactic 
bulge. 

The first measurements of microlensing were reported by Udalski et al. (1993) and 
Alcock et al. (1993). Several microlensing searches are currently underway, including 
the OGLE collaboration (Udalski et al. 1993), the MACHO collaboration (Alcock et al. 
1993), the EROS collaboration (Aubourg et al. 1993), the AGAPE collaboration (Alard 
et al. 1995), and a collaboration monitoring M31 (Crotts et al. 1992). Groups that track 
microlensing events, once announced, include the PLANET collaboration (Albrow et al. 
1995), and the GMAN collaboration (Becker et al. 1997). 

Combined, microlensing search collaborations currently monitor tens of millions of 
stars down to a visual magnitude of about 19.5 over a small fraction of the sky, currently 
reporting over 100 microlensing events (for a good review, see Paczynski 1996). Of the 
first 45 gravitational lens candidate events toward the Galactic bulge, the brightest star 
recorded to undergo a microlensing event by the MACHO collaboration was of unlensed 
visual magnitude 16.4 (Alcock et al. 1997). 

A primary motivation for this study was curiosity about how bright a star was 
undergoing a detectable microlensing event at any given time. Previously, Colley & 
Gott (1995) estimated that a microlensing events becomes bright enough for naked eye 
observation nearly every 5000 years. Perhaps, however, a star that is significantly brighter 
than those being monitored is currently undergoing detectable microlensing. Microlensing 
of this star might be relatively easy to track, as it would be relatively bright and in a 
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less crowded field. This paper is an attempt to estimate the magnitudes of the brightest 
microlensing events on the entire sky, and to discuss the possible advantages of finding 
these events. §2 will detail how this estimate will be made. In §3, results will be presented, 
while some discussion and conclusions will be given in §4. 

2. All-Sky Microlensing Rate Estimates 

Nemiroff (1986) gave a first estimate of the all-sky microlensing rate. There it was 
shown that to amplify source light at least A absi a lens with Schwarzschild radius Rs must 
fall into a "detection volume" V = (2/3)nRsDQ S § abs , where distance D has subscripts O, 
L, and S which stand for observer, lens, and source, respectively. & a bs = A/\/A 2 — 1 — 1, 
and here carries the subscript "abs" which refers to absolute amount of source magnification. 

Assuming all stars are identical and uniformly distributed, the likely number of lens 
events N above absolute amplitude A abs out to distance D os is found by summing the 
space internal to all the lens detection volumes, so that N = (8/15)7r 2 n 2 Q a b s RsDQ S , where 
n is the number density of stars. This early estimate did not include the brightness of the 
lensing star itself. Estimates of the effects of the brightness of the lens or an unlensed 
component have been done previously with regard to existing microlensing surveys by 
Kamionkowski (1995), Nemiroff (1997) and Han (1998). Nemiroff (1997) has shown that for 
a field of lenses with identical brightness, the apparent optical depth is a factor of about six 
less than the absolute optical depth, which includes the assumption that the lens is dark. 
Including this factor in the previous estimate, the number of lenses expected above the 
apparent amplification factor of A app would be 

N ~ (A/Ahyn 2 § app R s Dl s . (1) 
Given n = 0.1 stars pc~ 3 and R s = 3 km, <& app = 1/2, the single closest lensing 
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event is expected at D min = 1200 pc. The total number of stars in this volume is 
iV(stars) = (4/3)irDQ S n = 7.0 x 10 s . Given the canonical star is like the Sun and has 
absolute visual magnitude My ~ +5, its apparent magnitude at 1200 pc is an estimate of 
the brightest microlensing event occurring on the sky at any one time: my ~ 15.4. 

Although instructive in calibrating intuition, this estimate did not include the 
anisotropic and inhomogeneous stellar distributions that surround our Earth and our 
Galaxy. A more accurate estimate might be significantly different. Toward this end, an 
all-sky microlensing rate estimate is made here by employing the Galactic model published 
originally by Bahcall & Soneira (1980) and subsequently updated by Bahcall (1986). The 
Bahcall-Soneira Galactic model obtains a reasonable fit to the observed star counts with 
a relatively simple disk and spheroid model. Bahcall (1986) demonstrates accuracy of the 
model in several directions down past 20th magnitude, although all of these directions are 
at least 20 degrees away from the Galactic plane. The model also included large-scale 
affects of reddening and obscuration by Galactic dust. 

A FORTRAN-program version of the Bahcall-Soneira model was downloaded from 
J. Bahcall's web site in December 1997 and adapted to address microlensing. The 
Bahcall-Soneira code tracks both main-sequence and giant stars for both a Galactic disk and 
spheroid. For main-sequence disk stars, the code was run with absolute visual magnitude 
My cut-offs of —6 at the bright end at +16.5 at the dim end. For disk giants, My was taken 
to be —1.5 at the bright end. For spheroid stars, My was taken to be —3 at the bright end. 

In general, the total number of lensing events visible is given by 

Mvr r D oi X n rM BRT rD OS [M BRT 

N= dn D os dD os / n(M) dM s / dD OL / n(M) irb 2 L dM L , (2) 

Jo Jo Jm dim Jo Jm D im 

where Q is a solid angle on the sky, Dq^ x is the maximum source distance, Mj^m is the 
cutoff for absolute visual magnitude at the dim end, Mbrt is the cutoff for My at the 
bright end, M s is the absolute visual magnitude of a source star, and M L is the absolute 
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visual magnitude for the lens. The impact parameter bi is the furthest distance from the 
observer-source axis a lens could be placed, at distance D OL , and still amplify source light 
by the amount A app needed for detection. From Nemiroff (1986, 1989, 1997) this impact 
parameter was found to be 



The brightness of the stellar lens itself works to dampen the magnitude of observed 
microlensing events. More precisely, to achieve an apparent lensing amplification of A app , 
the lens must create an absolute amplification A a b s > A app to be seen over its own brightness 
(Nemiroff 1997): A abs = A app {l + Il/Is) ~~ Il/Is, where / indicates apparent brightness. 

Microlensing rates from several candidate scenarios were made. In the first, it was 
assumed that all stars brighter than a given unlensed magnitude limit would be monitored 
for microlensing, a case referred to as "without magnification bias." In the second scenario, 
it was assumed that a star more faint than a given magnitude limit could be amplified 
above the limit and be detected. The second scenario therefore includes magnification bias 
(Nemiroff 1994). In the third scenario, only stars in the Galactic disk of the Bahcall-Soneira 
model were included, without magnification bias, to give a indication of the angular 
distribution of the microlensing events. 

In all of the computations, giant stars were not allowed to be lenses, since one or 
both amplified images might impact the lens surface and be absorbed (Ftaclas 1998). 
The mass of main sequence stellar lenses was taken from Mihalas & Binney (1981) to be 



Equation 2 above was used to compute N(m), the number of expected lens events 
anywhere on the sky. For any magnitude limit m, however, the maximum expected 
amplification Auax on the sky can be computed directly from N{m) since N{m) scales as 
$ app . In fact, when Q app = 1/2, N — r, the optical depth (Vietri & Ostriker 1983; Nemiroff 




(3) 



M/M Q = (L/L e ) 



0.3125 
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1989). When N events are expected at the level of $ app = 1/2 (where A app = 1.34), then 
fewer events are expected at a higher A app (Nemiroff 1997). Specifically, the level of & app 
where one expects a single event is at $ opp /(l/2) = N. The corresponding amplitude is 

' N + 2 (4) 
VN 2 + AN y ' 



3. Results 

Figure 1 shows a plot of the expected optical depth to stellar microlensing as a function 
of apparent visual magnitude my. When the inverse of the optical depth becomes on the 
order of the number of stars being observed, then there is a 63 % chance that at least one 
of the stars is amplified by a factor greater than 1.34, and detection of gravitational lensing 
becomes likely. Including stars in both the Galactic disk and the central spheroid, stellar 
microlensing becomes likely between 16th and 17th magnitude. The Bahcall-Soneira model 
used predicts almost exactly 100 million stars exist over the entire sky that are brighter 
than 17th magnitude. 

Perhaps unexpected is that Figure 1 shows that optical depth can actually be a 
decreasing function of the survey apparent magnitude limit. This is caused by the shape 
of the stellar luminosity function, which may cause a disproportionate number of stars 
in a given magnitude range to be relatively close to the observer, hence contributing 
relatively little to the optical depth. This decrease in optical depth means only that the 
probability of detectable lensing per star has decreased. Since star counts increase greatly 
as magnitude dims, the total probability of finding microlensing somewhere on the sky at 
fainter magnitudes still sharply increases. 

Figure 2 shows plots of the expected maximum amplitude of the brightest lens 
event on the sky as a function of the faintest apparent visual magnitude monitored. 
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Inspection of this plot again indicates at about 17th magnitude, an A Max > 1.34 event is 
expected, even without magnification bias. A search sensitive to magnification bias could 
discover that a star temporarily brighter than 16th magnitude is actually being detectably 
microlensed, somewhere on the sky. The more stable the relative photometry, the smaller 
the microlensing enhancements it would be sensitive to, the brighter the star that would 
show detectable microlensing. Alternatively, were all stars brighter than 20th magnitude 
monitored, at least one would be expected to be undergoing an A > 3.2 event at any given 
time. 

Figure 3 shows the number of A > 1.34 lensing events expected over the sky as a 
function of visual magnitude. This plot once again indicates a 63 % chance that one star 
with unlensed magnitude brighter than 17 is being amplified by A > 1.34, and that it is 
likely a star appearing to be brighter than 16th magnitude is undergoing an A > 1.34 
microlensing event. 

Additionally, Figure 3 indicates that on average, two stars with unlensed magnitude 
brighter than 18 are undergoing detectable microlensing, while potentially 9 stars brighter 
than 20th magnitude are expected to undergo simultaneously detectable gravitational lens 
magnifications with A > 1.34. 

Figure 4 shows a contour plot of the number of stars per square degree brighter than 
17th magnitude expected to show detectable microlensing. Figure 4 is plotted in Galactic 
coordinates and is centered on the Galactic Center. The plot shows that the most probable 
place to find gravitational lensing in our Galaxy is first toward the Galactic center, and next 
in our Galactic plane. In the direction of the Galactic center, the success rate of a lensing 
search monitoring all stars brighter than 17th magnitude would be on order 3 x 10~ 4 lenses 
degree -2 . 

These results estimate the static probability of gravitational lensing - assuming stars 
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don't move. But stars do move and lensing events have finite durations. The methodology 
for making dynamic lensing probability estimate, allowing for relative lens, observer, and 
source motion, was discussed previously by Nemiroff (1991). 

Taking the detection volume defined above and allowing it to move will allow it to 
"catch" increasingly more lenses as a function of time. Were the source to move with 
transverse speed v relative to the observer for a time t, it would sweep out a volume 

dynamic = (l/2)lTabvt = (tt/4) R^ D^^Vt, (5) 

were a = D Q s/2 is the semi- major axis, and b = yj R s D s^ap P is the semi- minor axis of the 
ellipsoidal detection volume. 

Given canonical values of R$ = 3 km, D os = 1 kpc, $ app = 1/2, v = 100 km sec -1 , 
and t — 3.1 x 10 7 sec, the ratio of dynamic to static volumes is about 5. This means that 
the probability of lensing is about 5 times as great over the course of a year than during a 
single observation. 

In general, the duration of a microlensing events would be on order tdur ~ 2b/ v. Given 
the canonical values listed above, a typical duration for a microlensing event would be 
tdur ~ 50 days. The brightness of the lens and the maximum amplification of the event is 
automatically included in § ap p- In general, higher amplification events, closer sources, and 
brighter lenses would cause shorter duration events. 

From Figure 3, the magnitude at which N = 1/5 is about 15th magnitude without 
magnification bias, and about 14th magnitude with magnification bias. Therefore, were all 
stars brighter than 15th magnitude monitored for a year with observations spaced shorter 
then the event duration, at least one is likely to show an A > 1.34 event. 



4. Discussion and Conclusions 



Directly applying stellar densities derived from the Bahcall-Soneira model to 
microlensing leads to several inaccuracies. For example, this model is not designed to 
directly indicate the space density of binary stars, dark matter lenses, or the Galactic bar, 
each of which could alter the microlensing rate. No corrections for these were made in the 
present calculation, although it is believed that these factors would only increase the actual 
observed lensing rate. 

One conclusion that can be drawn is that microlensing effects are present in any set 
of 100 million source stars that can possibly be monitored. In fact, the brightest 100 
million stars are among the least likely set to show microlensing, since their proximity 
allows relatively little room to place the lenses. Additionally, these stars are so close 
that wide-angle measurements are needed just to monitor them. Searching for all sky 
microlensing, however, should yield the closest lenses and the brightest sources. 

The relative brightness of these microlensing events might make them more than a 
curiosity. At 15th magnitude, most observatories would be capable of follow-up coverage of 
an announced, on-going microlensing event. This could lead to microlensing events with 
extraordinarily good time coverage. Good time coverage, in turn, creates an increased 
chance of recovering brief bumps in microlensing light curves caused by planets orbiting the 
stellar lens (see, for example, Mao & Paczynski 1991; Gould & Loeb 1992). 

The relative proximity of these lenses and/or sources might be expected to create an 
increased chance that parallax could be measured for them, and hence foster the ability to 
better understand lens and source distances, masses, and even transverse velocities. For 
lenses crossing the disk of the source, the lens will act like a huge magnifying glass being 
swept across the source star face. Potentially recoverable information includes the radius, 
shape, and brightness distribution of the source star (Nemiroff & Wickramasinghe 1994; 
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Gould 1994; Witt & Mao 1994; Gaudi & Gould 1998). 

The relatively high angular separation of these microlensing events might be expected to 
allow optical depth estimates toward widely separated directions in our Galaxy. Therefore, 
event rates would give information on the Galactic stellar and dark matter distributions. 
The optical depth to lensing has yet to be accurately predicted in any direction - the result 
is always a surprise. Possibly optical depth estimates along new directions would also yield 
new surprises. 

Some might consider the above results discouraging, indicating that the number of stars 
is too high, the event rate too low, the data rate too high, and the monitoring instruments 
needed are too unusual to search for microlensing over the entire sky. Others might consider 
that the technology need to carry out a large area search capable of detecting microlensing 
is not only possible, but surveys that could be considered prototypical are already being 
implemented. Examples of such monitoring systems include LOTIS (Park et al. 1997) 
and ROTSE (Marshall et al. 1997), originally created to search for optical counterparts to 
gamma-ray bursts. LOTIS is already generating daily images of approximately n steradians 
down to 13th magnitude while searching for the optical counterparts of gamma-ray bursts. 
A planned improvement creating "SUPER-LOTIS" will scan the sky every three weeks 
down to a limiting visual magnitude of 19, at which time bright star microlensing events 
would be expected to exist in its database. 

Lastly, it should be pointed out that the "secondary" information on the intrinsic 
variability of bright star systems gathered in a wide-angle microlensing search should be 
of significant value. Interesting causes of variability includes stellar pulsations, binary star 
eclipses, and recurrent novae. Arguably, since transient events of even unknown cause could 
be followed by small telescopes common both to professionals and amateurs, "triggers" 
created by a wide-angle monitor could be of interest regardless of whether they were related 
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to microlensing. 
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Fig. 1. — A plot of the optical depth to gravitational microlensing averaged over the entire 
sky, versus apparent visual magnitude. The solid line depicts the optical depth to a model 
incorporating both the Galactic disk and the central spheroid, while the dashed line depicts 
optical depth including only the disk. 

Fig. 2. — A plot of the maximum expected gravitational-lens induced amplification expected 
over the entire sky, as a function of the faintest visual magnitude observable. The solid line 
depicts a survey of a fixed sample of stars so that no magnification bias should operate. 
The dashed line depicts a survey where stars below the survey magnitude limit could be 
micro lensed above the survey limit: magnification bias. The dot-dashed line depicts a survey 
of only disk stars without magnification bias. 

Fig. 3. — A plot of the number of lensing events with apparent amplitude A > 1.34 expected 
over the entire sky at any time, as a function of the faintest visual magnitude observable. 

Fig. 4. — A plot of our Galaxy contoured by the expected surface density of microlensing, in 
number of events with A > 1.34 per square degree, without magnification bias. The Galactic 
center is at the origin. 
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